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' Abstract 

-xt- : 

Resummation in QCD provides insight into the evolution of final state jets from short 
to long distances, and of accompanying interjet radiation. Applications to event shapes, 
including the recently-proposed angularities, suggest experimental tests of the interrelations 
between weak- and strong-coupling dynamics. 
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: 1 Introduction 

> ■ 

^ ■ Jets are the imprint on final states of dynamics at short distances, whether from transitions 
within the standard model (Drell-Yan annihilation or QCD scattering, for example) or from 
the creation of new degrees of freedom. An individual jet, however, provides scant clues to 
its origin. Indeed, a jet is simply a collection of co-moving particles in relative isolation. The 
number of particles in a given jet, and consequently its total energy and momentum, are never 
uniquely defined. Nevertheless, the distributions of particles and momenta within and between 
high energy jets encode information at all scales, from the largest momentum transfer or mass, 
through QCD evolution to the scale of the strong coupling, Aq CD . From the point of view of 
quantum field theory, the system passes through a stage of weak coupling at the shortest distances 
on to strong coupling and nonperturbative dynamics at the longest. Starting with perturbation 
theory, resummations provide one bridge between these regimes. 
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2 Resummation: Why? When? How? 



Only the most inclusive observables depend on a single hard scale. More common - and more 
interesting - are functions of two ordered but perturbative scales, Q ^> Q\ ^> A QCD . Varying the 
lower scale, Q\ allows us to move continuously between the ultraviolet and the infrared. Typically, 
as Qi decreases, the perturbative series develops one or two logarithmic enhancements in the ratio 
Q/Qi for every power of a s . These enhancements may often be organized, that is resummed, 
to all orders in perturbation theory. Examples include enhancements that are explicit in the 
perturbative cross section, such as those in the transverse momentum distributions (Qi = Qt) of 
vector bosons [1], and enhancements that are implicit in inclusive distributions, as in threshold 
resummation for inclusive Higgs production [2], integrated over Q T . 

Perturbative calculations in QCD depend on the fundamental property of asymptotic freedom. 
The most tractable perturbative quantities are one-scale cross sections that are fully infrared safe, 
that is, cross sections that can be expressed as series in a s with finite coefficients. Infrared safe 
observables are essentially descriptions of the flow of energy [3] . The class of infrared safe, single- 
scale cross sections is limited to e + e~ total and jet cross sections, however, and the latter are 
single-scale only when the masses of the jets are comparable to the total energy Q. When jet 
masses become small compared to Q, it is necessary to resum, even in these infrared safe cross 
sections. 

Resummation of two-scale logarithms can be derived whenever a cross section (sometimes 
amplitude) is a product or convolution of factors that separate the disparate scales Q and Qi 
through the introduction of a third scale, the factorization scale, \i ^> Aq C d- Schematically, a 
factorizable cross section a that depends on a very large momentum scale Q and a much softer 
scale m can be written as <j(Q, m) = uj(Q, //)<£)/(//, , m). Here m labels the soft scale(s), typically 
light quark masses and Aqcd 5 although it may also represent the lower of two perturbative scales 
(say a jet mass) in an infrared safe cross section. Whenever there is such a factorization, there is 
evolution. Since the physical cross section cannot depend on the factorization scale, the variation 
of the short distance function, u, with /i must cancel that of the long distance function, /, 

d n x dlnf dlnu 

/i— lno- phys (Q,m) = =>> = - P{a 8 W) = > t 1 ) 

where the "kernel" P can depend only on the variables that the two functions hold in common. 
Wherever there is evolution there is resummation [4], which is simply the solution to the evolution 
equation or equations. An alternative route to resummed cross sections is based on coherent 
branching, which analyzes repeated gluon emissions [5]. One may think of each branching as a 
"minifactorization" , an incremental step from short to long distances. 

Factorization proofs [6] are based on the quantum incoherence of dynamics at incommensurate 
length scales, and of the evolution of systems of particles mutually receding at nearly the speed 
of light. These considerations are reflected in the structure of factorizable cross sections near an 
"elastic" limit, with final states that may be characterized by a single hard scattering and a fixed 
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number of jets, 

Wjets 

a(Q, a + b^ iV jcts ) = H ® TV/a ® 7V/ 6 <E> II J * ® 5 ■ ( 2 ) 

i=i 

The convolutions may be in partonic momentum fractions, transverse momenta or energies, 
depending on the observable. Physics at the hard scale is in H, while the remaining functions 
generate perturbative logarithms and nonperturbative dependence. In effect, the calculation of 
any such cross section can be broken down into a set of standard components: the product of P's, 
which describe how two partons a' and b', which result from the quantum evolution of partons a 
and b, collide at H, leaving behind forward jets of "spectators" , and producing a set of outgoing 
jets, Ji and coherent soft emission S. For example, in hadronic collisions when partonic energy 
is just sufficient for W or Z production, the cross section takes the form H <g> V q / a <E> Vq/b <S> S. 
Such an inclusive cross section is sensitive to long distance dynamics only through the forward 
jets and the soft radiation. Similarly, for e + e~ — > 2J the cross section can be factorized into a 
hard part times two jets, H <g) J q ® Jg <g> S ', while a DIS structure function near x — 1 breaks 
up into H ® V q /a ® J q ® S. Because the cross sections factorize in these limits, their logarithmic 
dependence on the masses of the jets and on 1 — x, respectively, may be resummed [4, 5, 7]. 



3 Application: Angularities in e + e 

Among the many applications of resummation, event shapes near the two-jet limit in e + e~ 
annihilation have received perhaps the most attention, in large part because of the large lever 
arm provided by the LEP data, both in the large scale, the total cm. energy, and in the smaller 
scale, typically the jet mass. 

A generalized class of event shapes, the "angularities" , were proposed in Ref. [8] (reanalyzed 
and renamed by Berger and Magnea in [9]), with the motivation of providing a new parameter 
that interpolates between distinct traditional measures of jet substructure. For final state N, we 
define 

T - = 7j E Msm^ra-lcos^l) 1 -"^ EiW^cosOi). (3) 

^ iinN y iinN 

Here $i is the angle between the direction of final-state particle % and the thrust axis. Special 
cases are a = 0, the thrust, a = 1 broadening, and the total cross section, a — > — oo. The "elastic" 
limit for any finite a is r a = 0, where the final state consists of two perfectly collimated jets. 
Following the reasoning outlined above to next-to-leading logarithm (NLL) in r a , the differential 
cross section factorizes in the space of Laplace moments [8], 

a(T a ,Q,a) = cr tot / due VTa [Ji(v,pji)} 2 , 
J c 

Jfapji) = [dr a e-^ J i (Tj i ,p Ji )=e^ E ^> a \ (4) 
J o 
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into a product of jet functions, J. At the level of NLL, it is possible to define S cS = 1, which 
essentially serves as a definition of the jets. This definition is equivalent to the calculation of the 
cross section in terms of jets evolving independently according to coherent branching [5]. 

Logarithms in the transform variable v are in close correspondence with those in r a . In 
transform space when a < 1, logs of v exponentiate in Sudakov form, with up to n + 1 logarithms 
at order a™ in the exponent E in Eq. (4), given in terms of known anomalous dimensions A(a s ) 
and B(a s ) by 

1 , ,- uQ 2 

du r — 



E(v,Q) = 2/^ / d J^A{a s {p T ))U-^T/QY _ A 



+ 



l -B{a s {^Q)) ( e -^ 2/(2 - a) - 1 



(5) 



A characteristic feature of resummed cross sections, illustrated by this expression, is an integral 
over scales of the running coupling. Taken literally, these expressions are ill-defined, from regions 
where the integration variable px is of order Aq^. Although this singularity can be avoided 
while retaining NLL accuracy [7], it is also useful to study the implications of such ambiguities, 
as we shall see below. 

The a-dependent expression (5) has yet to be confronted with real LEP data, except for the 
case a = 0, the thrust. Fig. 1 shows data for the closely-related heavy jet mass distribution at 
the Z pole [10]. The perturbative-only NLL prediction has the right shape overall, but is shifted 
toward smaller values of jet mass. We will attribute this shift to nonperturbative corrections 
below. 




Figure 1: Heavy jet mass distribution at the Z [10]. Dashed line: NLL resummed. 



These considerations may be generalized to jet shapes in deeply inelastic lepton-hadron and 
hadron-hadron scattering when the overall final state can be characterized by a definite set of jets 
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accompanied by soft radiation, as in Eq. (2) above. Such cross sections, were termed "global" 
by Dasgupta and Salam in Ref. [11]. Recently, Banfi, Salam and Zanderighi have extended NLL 
resummation to a wide class of global observables in e + e~, DIS and hadron-hadron scattering, 
by developing an innovative software package [12]. 

4 Non-global Logs: Color and Energy Flow 

Complementary to jet shapes are descriptions of interjet energy flow. A simple illustration is 
shown in Fig. 2, where we trigger on two jet events in the scattering of particles A and B, and 
measure the inclusive distribution En(£'), where E > En > 0, with En the energy that flows 
into some angular region Q, away from both the collision axis and the jet directions. Quantities 
like Hq{E) are sometimes referred to as radiators. 




Figure 2: Geometry for energy flow observables. 

We can imagine (at least) two choices for such a cross section. First, it may be fully inclusive 
in the region Q between Q and the jets. In this case, the number of jets is not fixed, and the 
observable is nonglobal in the terminology Dasgupta and Salam [11]. This observable cannot be 
factorized into a fixed number of jets as in Eq. (2), and as such cannot be resummed to a simple 
exponential in the same way as the event shapes described above. Alternatively, we may limit 
radiation into region by constructing a correlation with an event shape such as r a that fixes 
the number of jets [8, 13]. 

Cross sections where the number of jets is not fixed are not fully understood, but they remain 
infrared safe, so that we should be able to learn about them in perturbation theory. Indeed, 
Banfi, Marchesini and Smye [14] showed that at leading logarithm, a™ \n n (\/~S / En) , and in the 
limit of large numbers of colors, N c , these cross sections obey a beautiful nonlinear evolution 
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equation, 

d A Z ab (E) = -d A R ab H ab {E) + / dN ab ^ k {H ak {E)ll kb {E) - Kb(E)) , (6) 

Jk not in Q 

where d& = E8e- Here dN ab ^ k describes the angular distibution of radiation from a pair (dipole) 
of color charges, and R ab the corresponding logarithmic angular integration, where the (3 are four- 
velocities of partons within the jets that radiate soft gluons directly into VL (case 1 in Fig. 2), 

dn k (3 a -(3 b f QdE> r 

dN ab _ k = — Tr-Tnr-T, R <* = L-et L dN ^ ■ ( ? ) 



4tt Pk-PbPk-Pa Je E< Jn 

These quantities control the linear term on the right-hand side of Eq. (6). The nonlinear terms 
describe radiation from a hard gluon of momentum k into region Q. This gluons acts as a new, 
recoil-less source of further emission into region Q (2 in Fig. 2). In the large- N c limit, any gluon 
may be thought of as a pair of sources q(k)q(k), with quark- ant iquark color. In this way, for 
large N c , the combination q(a)G{k)q{b) is equivalent to an independent pair of dipole sources, 
q(a)q{k) © q{k)q(a). An intriguing relation has been pointed out of this equation with small-x 
evolution in the unitarity (saturation) limit [15]. 

In the more restrictive approach, the correlation with an event shape, for example an angu- 
larity, can fix the number of jets by setting r a Q ~ E [8], and we may factorize and resum as 
above, da/dEdr a ~ S(E/r a Q) <g> da vcsuu J 'dr a . In a combination of both approaches [13], one may 
trigger on two narrow jets, and take the limit E/(r a Q) — > 0, using the nonlinear evolution above 
for the function S that describes radiation into the interjet region. 

With methods such as these, we can study the influence of color flow at short distances on 
energy flow at wide angles [16], including applications to final state rapidity gaps [17]. 

5 From Resummed Perturbation Theory to 
Nonperturbative QCD 

We now return to the interpretation of expressions like (5), where the argument of a s vanishes in 
the infrared limit, As required by infrared safety, if we reexpand the running coupling in terms of 
a s (Q) for any fixed momentum scale Q, the result is finite at all orders. The divergence associated 
with the running coupling comes from an n\ behavior of this expansion at order a n s [18]. We 
treat such contributions as ambiguities in the perturbative expansion that may be resolved by 
supplementing the series with nonperturbative parameters, which turn out to be associated with 
power corrections in the large scale Q. Shape functions [19] organize the dominant corrections 
for event shapes such as angularities to all powers of Q. 

To be specific, in Eq. (5), we treat the region > k using perturbation theory, with k > Aqcd 
a new factorization scale. For p T < k, we expand the integrand and replace the complete series 
of powers of wp T jQ by J^np, the shape function, 

n k 2 




E(u,Q,a) = E PT+T ^- £ ^(-^1 f^p n T A(a s (p T ))^E PT + \nf^J^, K )(8) 



o 
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The jet function is now factorized in z/-space into perturbative and nonperturbative functions, 
resulting in a convolution after the inverse transform (4), 



Q) = r Q UMO °Mr a Q - £, Q) ■ 
Jo 



(9) 



Because the momentum space shape function is independent of Q, once it is chosen to describe 
the r a -distribution at some fixed Q, say the Z mass, it gives predictions for all Q [9, 19, 20, 21]. 

Fig. 3 illustrates for the heavy jet mass the quality of fit that may be achieved in this way for 
a very wide range in Q. What we learn from such an event shape is illustrated by the functional 




Figure 3: Jet mass fit with a shape function at the Z mass, with predictions and data for other 
energies. 



form discussed in [20] /o,np(p) = const p e~ p . The parameter c may be interpreted in terms 
of the transverse momenta emitted per unity rapidity range, while d is related to the flow of 
energy between the hemispheres associated with the jets. More generally, shape functions are 
related to correlations of energy flow. This connection may be made in a manner that recalls the 
moment analysis of multiplicity distributions [22]. We introduce an energy flow operator £(Q) 
by its action on states, 



N 



smh, ...k N ) = J2 ^(Q - n^ih, ...k N ). 



(10) 



3=1 



For r a — > 0, the final state is characterized by two narrow jets accompanied by soft radiation. In 
this limit, moments of the shape function may be represented as [19] 



„ n 

dZCfaiO = I II <Mi w a (cose j )(S(Q 1 ) . ..S(Q n )) 



'IF 



where the angular functions w a (cos9), have been defined above in Eq. (3). The expectations in 
(11) are taken in the presence recoilless color sources (Wilson lines) that accurately represent the 
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coupling of the soft radiation to the jets. The logarithm of the shape functions is a cumulant 
expansion in these expectations [20]. 

By generalizing the thrust and related event shapes to angularities, we recognize in Eq. (8) an 
interesting scaling property for the associated event shapes event shapes, which may be thought 
of as a test of the rapidity-independence of nonperturbative dynamics [9, 21], 

i 

— a 

(12) 

It would be interesting to confront this prediciton with LEP data. For the present, however, 
we must content ourselves with a comparison to PYTHIA, which is fairly encouraging, as Fig. 4 
shows [21]. It is worth emphasizing that most event shapes were invented for the e + e~ jet physics 






Figure 4: Scaling predictions for the shape function in a as tested by PYTHIA. 

of the late 70's, and there is more to learn by addressing existing data with new analyses [23]. 

Also relatively unexplored are power corrections for single particle inclusive cross sections. 
Here we would like to understand single-particle cross sections from fixed target to collider 
energies. Using the formalism of Ref . [24] , and the analysis of resummed cross sections described 
above, it is possible to study the energy dependence of power corrections in cross sections at 
measured xt = 2pt/\TS [25]. For example, in the phase space limit, xt — > 1 for A + B — > 7 + X, 
the cross section may be written as the inverse transform of the Mellin moments of the fixed-order 
cross section with a resummation of higher order logarithms, in a form analogous to (5), 

d(JAB- 

— .-^ 1 

-ioo 2ni 



Pt 



dp 



T 



-JV-1 



e E PT (N,p T ) e SE NP (N,p T ) 



(13) 



Isolating low scales as above, we derive a nonperturbative exponent with the N dependence 

1 , / , / 4p; 



5E NP (N) 



N pt 
const. — =- in — 
pi N 



5 Enp(xt) = const. 



p\ In 2 



4p2 



In pt In 




(14) 



where we have used the conjugate relation of the variables iV and \ax T . We find a complex 
behavior in S at fixed pt, associated with energy conservation [25]. 
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6 Hopeful Conclusion 



Resummations can bring perturbative QCD to the doorstep of nonperturbative field theory. 
This analysis is still evolving, and alternative treatments of the perturbative/nonperturbative 
transition are possible. Theoretical and experimental studies of the interplay of color and energy 
flow in hadronic scattering should aid these developments. Eventually we will learn to translate 
fully the language of partons into the language of hadrons. 
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